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We show that the class of quantum baker's maps defined by Schack and Caves have the proper 
classical limit provided the number of momentum bits approaches infinity. This is done by deriving 
a semi-classical approximation to the coherent-state propagator. 
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I. INTRODUCTION 

The introduction of 'toy' mappings which demonstrate essential features of nonlinear dynamics has led to many 
insights in the field of classical chaos. A well-known example is the so-called baker's transformation 0]. Interest in this 
mapping stems from its straightforward characterization in terms of a Bernoulli shift on binary sequences. It seems 
natural to consider a quantum counterpart to the baker's map for the investigation of quantum chaos. Unfortunately, 
there is no unique quantization procedure, and hence, we must embrace the possibility of different quantum maps 
limiting to the same classical baker's transformation. 

Balazs and Voros were first to conceive a quantum version of the baker's map. This was done with the help of 
the discrete quantum Fourier transform. Subsequently, improvements to the Balazs- Voros quantization were made by 
Saraceno || , an optical analogy was found Q , a canonical quantization was devised |^] , and quantum computing 
realizations have been proposed ^ A quantum baker's mapping on the sphere has also been defined ||. More 
recently, an entire class of quantum baker's maps was proposed by Schack and Caves using qubits [ {Tof . The Balazs- 
Voros quantization is but one member of this class. 

The classical limit of the Schack-Caves quantizations is the subject of this article. We explicitly derive a semi- 
classical approximation for the propagator in the coherent state basis. This enables us to give conditions upon which 
the Schack-Caves quantizations will behave as the classical baker's transformation in the limit H — > 0. We find that, 
provided the number of momentum qubits approaches infinity, the semi-classical propagator takes the form 



(b\B\ 



d 2 W 
dadt 



■ exp 



W{b*,a)/2h 



exp 



{\a\ 2 + \b\*)/4h 



where \a) and \b) are coherent states on the torus, and W{b* a) is a classical generating function. Similar propagators 
have been encountered before using spin coherent states [O, [l^] , but all may be thought of as variants of those derived 
long ago by Van Vleck |l3| and Gutzwiller |uj . Semi-classical propagators play an important role in the path-integral 
formulation of quantum mechanics pE| and the related theory of periodic orbit quantizatio n |]l6|| . The latter has been 
investigated thoroughly for the Balazs- Voros quantum baker's map |l7], [l^, |l^, 20 ; 2l], 22, 



In deriving a semi-classical approximation only for the one-step propagator, we avoid complications which will arise 
after many iterations of the mapping. For long time scales, simple quantum-to-classical correspondences will break 
down Q and one must incorporate the theories of decoherence |27| [2|, or continuous measurement ||(| . The 
classical limit of the Schack-Caves quantization has already been investigated p2| in this light using a decoherent 
histories approach |33|, [34|, [35|. However only a special case (9 = in our notationjwas considered. We proceed under 
an assumption that provided our one-step propagator agrees with the baker's transformation in the semi-classical 
limit h — > 0, decoherence will restore quantum-to-classical correspondences for long time scales. 

The paper is organized as follows. In Section II, we introduce the baker's map, both in classical and quantal form. 
Coherent states for a toroidal phase space are also introduced. In Section III our core results are presented. Here we 
derive semi-classical approximations to the coherent-state propagator and give conditions for when the Schack-Caves 
quantizations have the proper classical limit. Finally, in Section IV, we summarize our findings. 
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II. A CLASS OF QUANTUM BAKER'S MAPS 

The baker's map is a standard example in chaotic dynamics. It is a mapping of the unit square onto itself in the 
form 

q n+ i = 2q n - [2q n \ (1) 

Pn+1 = { P n+l2q n \)/2 (2) 

where q,p € [0, 1), [^J is the integer part of x, and n denotes the n-th iteration of the map. Geometrically, the map 
stretches the unit square by a factor of two in the q direction, squeezes by a factor of a half in the p direction, and 
then stacks the right half onto the left. 

The map's action may be rewritten in terms of the complex variable a = q + ip, 

A generating function for this mapping (up to an arbitrary constant) is 

W(b*,a) = — (3b* 2 + 8ab* - 3a 2 ) + - f 1 + - ) ( a + ib* - - ) I a + a*J (4) 
10 5 V 2/ \ 2/ 

assuming a + a* is non-integer. The classical baker's map may then be rederived via the relations 

dw dw , 

w = b ^ = a - (5) 

Interest in the baker's map is due mainly to the simplicity of its symbolic dynamics. If each point of the unit square 
is identified through its binary representation, q = ■ s\S2 . . . = SfcLi s k%~ k and p = • sqS-i . . . = X)^Lo s -fc2 -fe ~ 1 
(sj G {0, 1}), with a bi-infinite symbolic string 

s = . . .s_ 2 s_is • S1S2S3 • • ■ (6) 

then the action of the baker's map is to shift the position of the dot by one point to the right, 

s — > s' — . . . s- 2 s-is si • S2S3 (7) 

For a quantum mechanical version of the map, we work in the Z?-dimensional Hilbert space, Hd, spanned by either 
the position states \qj), with eigenvalues qj = (j + (3)/D, or momentum states \pk), with eigenvalues pk — [k + a)/D 
(j, k = . . .D — 1). The constants a, (3 € [0, 1) determine the periodicity of the space: |oj+d) = e~ 2ma \qj), \pr+d) = 
e 2m P\pk). Such double periodicity identifies Hd with a toroidal phase space. The vectors of each basis are orthonormal 
(ijllj') = (Pk\Pk') = Sk.k' arid the two bases are related via the finite Fourier transform 

(q j \F D \q k ) = (q j \p k ) = -j=ei«»>. 



For consistency of units, we must have 2nhD = 1. 

The first work on a quantum baker's map was done by Balazs and Voros Their expression for the map was 
given in the form 

where F D / 2 is the finite Fourier transform acting on half of the Hilbert space. Later Saraceno |}| improved certain 
symmetry characteristics of the map using anti-periodic boundary conditions (a = (3 = 1/2). Finally, taking again the 
anti-periodic Hilbert space, Schack and Caves |10| introduced a whole class of quantum baker's maps for dimensions 
D = 2 N . 

For these cases, we can model our space as the product of N qubits with a binary expansion association 

\qj) = \xi) <8> \x 2 ) ® ■■•<£> \x N ) x; G {0, 1} (9) 



10 10 1 

(a) 9 (b ) 9 (c) q 

FIG. 1: The Husimi function for each partially Fourier transformed state (^|) when N = 2, and, (a) n = 2, (b) n = 1 and (c) 
n = 0. 

where j has the binary expansion 

N 

j = x 1 ...x N -0 = J2x l 2 N - 1 and q j = i±ll^. (10) 



Next we rewrite the quantum Fourier transform as 

\ Pk ) = F D \q k ) = -±= J2 W)®\x2)®---®\x N )e 2 ^ x l 2N (11) 



/ 2 N 

X\,...,XN 



where y = y% . . .ypf ■ 1 = k + 1/2 and x = x% . . . xn ■ 1 = j + 1/2. 

The connection with the classical baker's map comes from its symbolic dynamics. In the quantum case, a string 
is created through the partial Fourier transform G n . It is an operator which Fourier transforms the N — n least 
significant qubits of a state 

G n (\xi)®---®\x n )®\ai)®---ig)\a N - n )^ = |xi)(g>- • -<g> \x n )® -i== \ x n+i)®' • -® \x N )e 27Tiax/2N ~" (12) 

V Xn+1,-..,XN 

where a and x are defined through the binary expansions a = a\. . . ajv-n ■ 1 and x = x n+ \ . . .xn • 1. In the limiting 
cases, we have Go = Fn and Gn = il- 

The analogy to the classical case is made clear through the definition 

\a N - n ■■■a 1 »x 1 ...x n ) =G n \\x x ) ® ■ ■ ■ (g) |x n ) (g) |ai) (& • • • ig) |ajv-n)J • (13) 

These states form an orthonormal basis and are localized in both position and momentum. They are strictly localized 
in a position region of width 1/2™ centered at • x\ . . . x n l, and are roughly localized in a momentum region of width 
l/2 N ~ n centered at 

Using this notation, Schack and Caves defined a whole class of quantum baker's maps B n (n = 1, . . . , N) 

B n = G n _ x o G^ 1 = ^ \a>N-n ■ ■ ■ aixx •x 2 ...x n ) (a N ^ n . . . a x • x x x 2 ...x n \. (14) 

X\ ,...,x n 

fll ,...,a,N—n 

The Balazs-Voros-Saraceno quantum baker's map is recovered when n = 1. In the language of Eq. @, we see that each 
quantum baker's map takes a state localized at lajv-n • ■ ■ a\»x\ . . . x n l to a state localized at lajv_ n . . . a\X\*x 2 ■ ■ ■ x n l. 
The decrease in the number of position bits and increase in momentum bits enforces a stretching and squeezing of 
phase space in a manner resembling the classical baker's map. In figures Ma), (b) and (c), we have plotted the 
Husimi function (defined below) for the partially Fourier transformed states (O) when N = 2, and n = 2, 1 and 0, 
respectively. The quantum baker's map is simply a one-to-one mapping of one basis to another. 
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It will be useful to rewrite our baker's map in the position basis. To do this, we first use ( fL2] ) and ( |l3| ) to rewrite 
Eq. (0) 



as 



" 2 N - n+1 



\*x%. ..x n zi . . .z N - n+ i)(»xi . . .x n yi . ..yN-n\ 



(15) 



Xl,...,x„ zl,...,ZN-n+l 
Oi,...,ajv-n yif-tVN-n 



where 



x exp 



N- 



2 N 



- ((j + 1/2) (I + 1/2) + 2 N - n Xl {l + 1/2) - 20' + 1/2) (k + 1/2)) 



AT-n+1 



j = a k 2 N - n ~ k , k = ]T y k 2 N - n - k and I = £ z fc 2 



AT-ra+l-fc 



fc=l 



fc=l 



fc=l 



Next, using (^), ( fl3| ) and the notation qj = (j + l/2)/D, q k = (k + X/2)/D, etc, we arrive at the quantum baker's 
map in the position basis 



B r . 



V2 



1 2"- n -l 2"-'* +1 -l 2"- 1 -! 



ON—n+l E E E E 

xi=0 j,k=0 



1=0 m=0 
N—n+1 _ 2" 



• |« + gm 2 Jv ^" +i - 2~ n )(qk + xi/2 + q m 2 
x exp [inD2 n (q jqi + 2~ n x 1 q l - 2qjq k ) 



N—n r>— n— 1 



(16) 



Note that it is possible to sum over the index j at this point. However the above representation proves to be most 
convenient when performing our semi-classical analysis. 
We will now introduce coherent states for Hd |3^, |37|, 38 : 



1 / 2 



TV \D 



1/4 oo D-l 



E E ex p 

(X=— oo j—0 

1 /2\ 1/4 ^ 



7rL> 



I a. 1 2 — a 2 ) — 7rZ? (oj — a + /i) 2 + Mr/i 



\Qj) 



?R (\af + a*) - nD (q* - 2 Qj 



9 [<£>(«&- a) |iD]|fc) 



(17) 
(18) 



where a = q + ip and #o is called a i/ie£a function J39| | 

oo 

o [z|t] = exp [i7r(r^ 2 + (2z + 



(19) 



The coherent states obey \a ± 1) = — exp[±7riDp]|a), |a ± i ) = — exp[=F7riI?g]|a), and are simply the standard (Weyl 
group) coherent states that have been (anti-) periodicized and then projected onto H.d- The normalization factor 
takes the form 

N 2 =e [qD\iD/2]6 [pD\iD/2] = 1 + 0(1/U) (D even) (20) 
and henceforth, will be set to unity. Finally, the Husimi function for our toroidal phase space is defined as | (-0 1 a,) | 2 . 

III. THE SEMI-CLASSICAL PROPAGATOR 

Our goal in this section is to explicitly calculate the semi-classical propagator in the coherent state basis. That is, 
we wish to obtain the leading term in an asymptotic expansion of the matrix element (b\B n \a) as D — > oo. Observe 
from Eq. flU) that in this limit, the total number of position and momentum bits N necessarily become infinite. 
However one has considerable freedom of choice on how this may occur (see Fig. ^). We wish to consider cases where 
the relative number of position and momentum bits approach infinity at different rates. To this end, we take the 
number of position bits to be in the explicit form n = n(N) = 8N + s, where < 8 < 1 is rational and s takes 



1 



2 3 4 5 6 7 

n (position bits) 



9 10 



FIG. 2: Different possible ways of taking the classical limit for the quantum baker's map. 



integer values. For ease of reading, we also introduce the constant <j) = 1 — 9 such that the number of momentum 
bits N — n = 4>N — s. We will also now identify the different quantum baker's maps through the new parameters, 
Bg tS = B n . The parameter 9 (<f) may be interpreted as the fraction of qubits allocated to the position (momentum) 
register as the total number of qubits N, is increased. In the analysis which follows we must consider the two cases 
9 = and 9 = 1 separately. The former contains the original Balazs-Voros-Saraceno quantization (n = s = 1) and 
will be investigated first. The second parameter s, describes an initial offset between the number of position and 
momentum qubits and has no semi-classical effect when 9 < 1. We will find, however, that s becomes important when 
8=1. 



A. Case 9 = 0. 



In this case the number of position bits remains constant n = s > 1 as we let D — ► oo. Using ( |f6| ) and ( p7| ) our 
matrix element becomes 

oo 1 D/S-l 2D/S-1 S/2-1 

(b\B , s \a) = sd-*i> E E E E E 

irD 



exp 



- — ( M 2 + H 2 ~ a 2 - b* 2 ) + m{n -v)-irD[q k + Xl /2 + (Dq m - 1/2)/ 'S -a + fi 



D[qi + 2(Dq m - 1/2) /S - b* + v\ + inSDlqjqt + x iqi /S - 2q 3 q k 



(21) 



where S = 2 s . To further the calculation, we now use variants of the Poisson summation formula to replace each sum 
with D in the upper limit, by an integral e.g. 



l/S l/S n/q , 

oo '„ y. oo JJ/b — 1 

E / exp[2ni(Dx-l/2)a]f(x)dx=- I E s{x-{j + l/2)/D)f{x)da=- E ( 22 ) 
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Re y 
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(a) 



(b) 



FIG. 3: Steepest descent paths for f(y). The original integration path along the real line (gray) is deformed to one where 
Im/(y) = constant (black). 



The result is 



1 5/2-1 1/5 2/5 

(b\B Q:S \a) = SD^ 2 J2 12 12 dx dy dz 

fj,,v,ot xi—0 m=0 

/3,7— — oo 

exp - ^- (\a\ 2 + \b\ 2 - a 2 - b* 2 ^j + in(^ -u-a-f3-j) 
- nD(y + xi/2 + m/S -a + fxj -itD^z + 2m/ S -b* + v} 
+ iirSD^xz + x\zj S ~ 2xyj + 2inD (xot + y/3 + z~f) 



(23) 



We are now ready to make a semi-classical approximation to our matrix element. More precisely, we will make 
a saddle-point approximation to the triple integral above. Only near a saddle-point will contributions from such 
an integral cancel the prefactor D 3 / 2 and lead to an 0(1) contribution for the matrix element. The saddle-point 
approximation can be written down immediately using well-known formulae found in any standard text ]40[ . However 
the limits in the above integrals are finite; therefore, the saddle point will not make a contribution in all cases. We 
need to consider this possibility carefully if we are to recover the classical baker's map. Hence, we will treat each 
one-dimensional integral separately and use the method of steepest descents. 
Consider first the y integration (with x a parameter) by defining 



1/5 

h= J dy exp [ - nDf(y)] 



(24) 



where 



f(y) = (y - A) 2 + 2i(Sx ~ 0)y 
A = a — x-l/2 — m/S — fx. 



(25) 
(26) 



An asymptotic calculation of this integral is enabled by deforming the integration path, currently along the real line, 
to one in the complex plane where Im/(j/) = constant. Two important cases are drawn in Fig. [|. The first (a) occurs 
when the saddle point (defined through f'(y') = 0) 



y = A - iSx + i/3 



(27) 
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satisfies < Key' < l/S. In this case, the steepest descent path is one which first travels along the hyperbola C\ 
from to — oo + i Im?/, then along the hyperbolic asymptote C2 to 00 + i Imy', and finally back to l/S via another 
hyperbola C3. Hence an asymptotic expansion for the integral I\ will be the sum of three parts, each associated with 
a different contour C\, C2 or C3 in the complex plane. Note that along the contours C\ and C3, the kernel attains 
its maximum at the end points and 1/5, respectively. Consequently, the leading term in an asymptotic expansion 
takes the form 



■1 



ttd/'( c ; 



exp [ - *Df(cj\ [1 + 0(1/Dj] 



(28) 



with c = or l/S. However the prefactor of 1/D in the above inhibits such terms from playing a role in the leading 
order approximation of our matrix element. As remarked before, we need prefactors of D~ x l 2 in each approximation 
of the three integrals in (^3|) in order to obtain an O(l) overall contribution for the matrix element. Hence we will 
simply discard the integration along contours C\ and C3, and make the approximation 

oo+* Im y' 

h » / dy exp [ - wDf(y)] if 0< Rcy' < l/S (29) 

— oc+i Im y' 

= D-^ 2 exp[~TTDf(y')]. (30) 

When Re y' < or Key' > l/S (Fig. ||(b)) the path of steepest descent no longer passes through the saddle point, 
and consequently, there will be no leading order contributions to the matrix element, i.e. we may set I± = 0. The 
third and final case occurs when Rey' = or l/S. One may investigate these possibilities by taking exactly one 
half of d30| ) as the approximation for I\. However, for simplicity, we will not deal with this case, except make the 
odd casual remark when needed. In summary, we take ( |30| ) as our approximation for I\ when < Key' < l/S, and 
otherwise zero. 

Similarly, for the z integration one has 



2/S 



with 



and the saddle point is 



I2 — J dz exp [ — irDg(z)j 


w D- 1/2 exp[-TrDg(z')] if 0<Rez'<2/S' 



g(z) = (z- B) 2 ~2i(Sx/2 + Xl /2 + j)z 
B = b* — 2m/ S - v 



z' = B + iSx/2 + ixi/2 + i-f. 

Now letting x vary again and setting 
h(x) = f(y') + g(z') - 2iax 

= ^S 2 x 2 + s(2iA ~iB + 2:1/2 - 2/3 + 7 - 2ia/S}x - 2iA(3 - iB{x x + 2 7 ) + [3 2 + (xi/2 + j) 2 

we have the final integral 

l/S 



JdxeM--Dh(x)] 



5S 2 D 



exp [ - nDh(x')] if < Kcx' < l/S 



with the saddle point at 



x' = -^(2iA-iB + Xl /2- 2/3 + 7 -2ia/S 



(31) 
(32) 



(33) 
(34) 



(35) 

(36) 
(37) 

(38) 
(39) 

(40) 
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Now, inserting our saddle-point approximations back into (|23j) and setting a = a\ + ia 2 and 6 = 61 + ib 2 , with a 
little algebra we obtain 



(b\B ,s\a) 



00 1 S/2-1 

r> E E E ex p 

fj,,u,a xi—0 rn— 
/3,7— — 00 



— ((2ox - &x - aix - 2/i + j/- 2a/S') 2 
5 I 



+ (a 2 - 26 2 + 33i + /3 + 2 7 ) 2 | + Mr(fz -i/-a-/3-7) + i7rD(aia 2 - &1&2) 

— |(2oi - 61 - x x - 2fi+ v - 2a/S){2a 2 + b 2 ) 

- (ai + 2h - Xx/2 - n - 2v)(xx + /3 + 2 7 ) + afo - 4/3 + 2 7 )/S'| 



provided that all three of the inequalities 



0< Rex' = — (^2a 2 + & 2 -xi/2 + 2/3-7j < - 
< Re y' = ~ fax + 26 x - x x /2 - fi - 2v) + Ua - 5m ) < ~ 



55' 







< \ Re z' = \ ( ai + 26i - xi/2 - /i - 2i/") - i f a + 5m) < i 
2 5V / 5i \ / i 



(41) 

(42) 
(43) 
(44) 



are satisfied. Otherwise the summand is taken to be zero. 

Note that although m no longer appears in the exponent, we cannot trivially evaluate the sum since not all values 
of m will satisfy ( fl3| ) and ((f4|). Consider the cases when the approximation ( f4l| ) becomes O(l). That is, 



2a l -b 1 -x 1 -2fi + v-2a/S = 
a 2 - 26 2 + x x + [3 + 27 = 0. 



Substituting ( fig ) into ( |42| ) we obtain 

< a 2 + (3 < 1 



or 



< 62-^1/2-7 < 1/2 



(45) 
(46) 



(47) 



and thus, the integers f3 and 7 give our periodicity in the momentum direction. Hence if we assume < a 2 ,6 2 < 1 
then we may set j3 = 7 = in (41), making note that we are discarding exponentially small Gaussian tails. Also note 
from (|47]) that we must have X\ = \2b 2 \ . 

Now, negating (Q) and adding it to (|4^) we immediately arrive at the inequality — 1 < a < 1. Hence we must set 
a = 0. This implies 



2m 2 ( , \ 

—g < g [ai + 2bi - xx/2 - fi - 2vj < 



2(m+l) 



from 



|), or equivalently 



< - ( ai + 26 x - xx/2 - n - 2v) < 1 



(48) 



(49) 



if we now drop the summation over m in (|4l|). Hence, following a similar procedure to the above, one can substitute 
( [lo] ) into the new inequality ( f49|) and deduce that under the assumption < a x , b\ < 1, the summand of (^l] ) becomes 
O(l) only when /i = t/ = 0. Furthermore, we will have xi = [2aiJ- 

The surviving term of the summation is our semi-classical approximation for the propagator: 



(b\B 0y 



■ exp 



nD 
~~5 



' C 9 9 

-|(2oi-6i- L2oiJ) +(a 2 -26 2 + L2aiJ) 
i(3aia 2 + 36i6 2 + 4aife 2 - 4a 2 6i) - 2i[2 ai \ (ax + 26i + 2a 2 + b 2 - [2aiJ/2) | 



•o(l) (50) 



where we have chosen xx — [2ax\ (and implicitly assumed ai 7^ 1/2 and b 2 =/= 1/2 by ignoring cases of equality in 
( fj^ - |44| )). All other terms in (4l|), being exponentially small, are discarded. 
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Note that the above approximation is 0(1) only when b is the iterate of a under the classical baker's map (j^). 
Furthermore, a little algebra reveals that our semi-classical propagator may be rewritten in the Van Vleck form 



(b\B Ja) = 



d 2 W 



dadb 



7 ex P 



irDW(b*,a) 



exp 



irD (\a\ 2 + \b\ 2 ) /2 + o(l) 



(51) 



where W(b* , a) is the classical generating function (Q). Hence we have shown that the class of quantum baker's map 
with 8 = will approach the classical baker's map in the limit D — > oo. 



B. Case < 9 < 1. 

We will now consider the case < 6 < 1. Using (|l^) and (|l7|) with n = 9N + s, our matrix element is 
S 



(b\Bg, a \a) = 



oo 1 D+/S-1 2D*/S-1 SD" / 2-1 

fi,v= — oo X!=0 j,k=0 2=0 m=0 



- ^- (|a| 2 + |6| 2 - a 2 - 6* 2 ) + nr(/x - v) - *£>((<# - l/2S)/D e + xi/2 + q 9 JS - a + ,i) 
((gf - 1/S)AD 9 + 2^/5 - 6* + i/) 2 + inSD*(q+qt + Xxqf/S ~ 2qfqf)'\ (52) 



exp 
- irD 



where again S = 2 s . Introducing the new summing variables q 6 m = q m D/D e , = qjD/D^, etc., enables us to convert 
the four finite sums over j, k, I and m, to integrals over x, y, z and t, respectively, using formulae similar to (p2 
The result is 



l/S l/S 2/S S/2 



(b\B e>s \a) = S^jDD^ jr ^2 J dx J dy J dz J 



dt 



exp 



7,«— — oo 



2 ^ - b* 2 



+ «7r(/i — ^ — a — /3 — 7 — k) 



- irDUy - 1/25) /L> e + x x /2 +t/S-a + ^j - nDUz - l/S) /D e + 2t/S - b* + 1/)' 
/5 - 2scy) + 2mD*{xa + y{3 + zj) + 2iirD e tK 

DD* °° 1 1 1 2 1/2 

I ^r— J dx J dy J dz J dt 



(53) 



exp 



li,y,a,fi xi=0 
7,k— — 00 

ttD 

~Y 

,2 



(|a| 2 + |6| 2 - a 2 - 6* 2 ) + wr(/x - i/-a-j9-7-«) 

- iiD(t + xi/2 - a + /i) 2 - nD(2t - b* + vf 

- 27rL> (y - 1/2) (t + xi/2 - a + fx) fS - 2ttL> (z - l) (2t - 6* + v) /S 
+ iirD^ (xz +xiz- 2xy) /S + 2inD' t ' (xa + y/3 + 27) /S 

- TrD*- e (y - l/2) 2 /S 2 - nD^ 9 (z - l) 2 /S 2 + 2mSD e tn 



(54) 



where we have rescaled the integration variables (x — > x/S, y — > y/S, z — > z/S and t — > then collected terms in 
the exponent with the same power of D. The terms with highest power are those containing t, and hence, we will 
consider the integration over this variable first. Define the integral 



1/2 

h= J dt exp -nD(t- A) 2 -nD(2t- B) 2 -2-kDHC + 2inSD e tn 




(55) 
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where the constants are 



A — a — xi/2 — jj, 
B = b* -v 

C = (y-l/2)/S + 2(z-l)/S. 



(56) 
(57) 
(58) 



We now wish to derive the contribution from I4 which gives the leading order approximation to our matrix element. 
This is done by taking a path of steepest descent for the function f(t) = (t — A) 2 + (2t — B) 2 . Note that when 0=1 
the third term of the exponent in (EjJ) also becomes dominant and must be incorporated into f(t). Hence the need 
to consider this case separately in the previous section. 

As before, we may discard all parts of our integration contour, except the segment (—00 + i Imi', 00 + i Ivat') which 
passes through the saddle point 



f 



A + 2B 



(59) 



It is only this contribution which will cancel the prefactor \f~DD^ in (|4|) to give an 0(1) overall contribution to the 
matrix element. Hence we make the approximation 



oc+i Im t' 



dt exp - nD(t - A) 2 - nD(2t - B) 2 - 2-kDHC + 2mSD e tK 



(60) 



— oo-H Im t' 

1 



'5D 



exp 



-^(2A- B) 2 - ———-(A + 2B)C + 2 -^(A + 2B)S. 
5 5 5 



-C 2 



ttD 



2 .,2 



-S Z K 



2in 



-CSk 



if < Ret' < 1/2, and otherwise zero. 

Substituting this approximation back into (|5~i| ) and simplifying we obtain 



(61) 



(b\Be tS \a) 



D4> 
755 



E E 

7,k— — 00 



112 

dx J dy J dz 




exp 



nD 



\b\ 2 -a 2 - b* 2 j + - v - a - - 7) - — (2A - B) 



-— (2A - B) (2y -z) + — — (xz +x lZ - 2xy) + — ^— (xa + y(3 + zj) 
-5-5- (22/ - z) - — (y + 2z) k + — ^- (A + 2B) - —^~S 2 k 2 



(62) 



The dominant terms in the exponent which contain the integration variables, are now those with as a prefactor. 
These terms do not define a saddle point, but instead, a line. Hence, it is advantageous to first decouple x, y and z 
in these terms using the following transformation 



(63) 



u 




" 1 


-2 


1 " 




X 


V 




1 


2 


-1 




y 


w 







1 


2 




z 



where the integration region (x, y, z) £ [0, ll x [0, 1] x [0, 2] is transformed to some parallel-piped f2. 
After making this transformation, Eq. (p2) may be rewritten in the form 



(b\B g , s \a) 



E E 



du dv dw exp 



"SW - one? ( M _ V > 



S 



20S 2 



F(w) (64) 
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(a) 



Re u 







u . 


(«,0) 


— > — ^ 








(«,0) 




/c 3 





(b) 



FIG. 4: Steepest descent paths for g(u). 



Re u 



where 



g(u) = ~-u 2 ~(E + ia)u 
h(v) = -v 2 + (E — ia)v 
E E 



\ (2A - B) + % - (xi/2 - 2/3 + 7) ee Ex + iE 2 

o 



(65) 
(66) 
(67) 



and 



F(w) 



exp 



itD 



(\a\ 2 + \b\ 2 -a 2 - 6* 2 ) + m{n - v - a - (3 - 7) - ^(2A - B) 



(A + 2B) Sk - ^—S 2 k 2 + |J (D* (x 1+ P + 2 7 ) - Sk) w 

Do V / 



(68) 



We have now arrived at a form where we can consider steepest descent paths for the integration variables u and 
v. Starting our program with the function g{u) and writing u in terms of its real and imaginary parts u = u\ + IU2, 
ones finds that the two hyperbolic asymptotes 



U\ — u 1 
1 



= u 2 -u 2 



U\ — u 1 = —112 + u 2 
are the steepest descent paths which pass through the saddle point 

U = u[ + iu' 2 = 2(iE — a). 



(69) 
(70) 



(71) 



However, only the first ( |69| ) can be used as an integration contour, since exp [— ttD^ Re g(u)/S] — > 00 on the other. 
Cases for when this asymptote (denoted by C2) is required to form part of the integration contour, and when it is 
not, are plotted in figures |^(a) and (b), respectively. Here the integration limits are denoted by the real numbers u 
and u, and need not be known explicitly for the moment. Note that C2 is included in the contour only when the 
intercept of the second asymptote (ww with the real line, denoted by u" , is between these limits. That is 



u < u 



2(Ex -E 2 -a) <u. 



(72) 



The importance of this inequality is not clear yet. We shall return to it after transforming back to our x, y and z 
variables. 
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The analysis of h(v) is very similar. In this case one finds that our steepest descent path will travel along the 
asymptote 



and hence, through the saddle point 



only when 



v\ — v[ = ~v 2 + v' 2 



■}[ + iv' 2 = 2(iE + a) 



v < v = v [ — v' 2 = 2(a — Ei — E 2 ) < v. 



(73) 
(74) 
(75) 



But what are our integration limits tt, u, v and v? Unfortunately, given the nature of the variable change, their 
values differ as one integrates over the volume element f2. Therefore, it is convenient to convert back to our original 
variables x, y and z which have independent limits. We know that whenever the point (u",v",w) belongs to our 
integration region 0, the two contours ( |69| ) and ( |73| ) will be included in our integration path of steepest descent. 
Therefore, by inverting our transformation (p3), we obtain 



u" = 2{Ei -E 2 -a) x" = -2E 2 

v" = 2(a-Ei-E 2 ) =► y" = |(a — E{) + -i 

w = w z" = — |(a— E x ) + 



and hence, the condition (u",v" 1 w) S 17 is equivalent to 



< x" < 1 
< y" < 1 
< z" < 2. 



Inequalities (|78|) and (79) may be rewritten as 



a — Ei < w < 5 + a — Ei 
-A(a-Ei)< w <5-i(a-E 1 ) 



(76) 



(77) 
(78) 
(79) 



(80) 
(81) 



which asserts that 

w<w<w (82) 

where w = max{a — Ei, — 4(a — Ei)} and w = 5 + min{a — E\, — 4(a — Ei)}. The requirement that w < w implies 

-l<a-Ei<l. (83) 

Having learned all we can from the method of steepest descents about the restrictions placed upon our integration 
parameters, we may proceed with the saddle-point approximation of (|64|). This is done by discarding all contours 
which do not make a contribution to the leading order approximation of our matrix element (e.g. for u, we discard 
C\ and C 3 in Fig. ||). Hence 



(b\B e ,s\a) 



£>4> 

exp 

1 /4 
5V 5 



u +ooe 



-r/4 



v + oce 



k/4 



E E 

7,re= — oo 



du 



dv / dw 



u — ooc 



t/4 



v — ooe 



9m ^-K v 



E E 

7,K — — oo 



S 20S 2 
dw exp 



F{w) 



(84) 



ttD 



|6| 2 - a 2 - fo* 2 ) 



nD 



+ m{li-v-a- fi-i) - — (2A-B) —Ea- 

S 



A-kD*~ 
55 2 



2mD" , A . „ irD e 
(A + 2B)Sk- — 



xi + (3 + 27) - S^jw 



(85) 
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Notice that we have actually made a saddle-point approximation about a line parametrized by w, over which, we still 
need to integrate. Currently there are no restrictions placed upon k. However the integral over w will be 0(D^^) 
unless k = D^(x 1 + (3 + 2 7 )/S. We may set k to this value by noting that we are only shifting the saddle point t' by 
the imaginary amount i"n{x\ + j3 + 2j) (see Eq. and thus, our approximation of I4 ( |6l| ) is unchanged. 

Hence, putting k = D^(x\ + (3 + 2j)/S and integrating over w, we arrive at 



(b\B 9 ,.\a) 



1 00 ^ 

eve E E 



I to — exp 



7 — — CO 



— I (2oi - 61 - .xi - 2/i + i/) 2 



+ (a 2 - 26 2 + xi + (3 + 2 7 ) 2 | + infyi -v-a-^-j) + i-KD(a x a 2 - &i6 2 ) 
-|(2ai - 61 - xi - 2/i+ v){2a 2 + b 2 ) 



2inD 



4ttD" 



- (01 + 2&i - xi/2 - /it - 2i/)(a;i + f3 + 2 7 )| 
provided that all three of the inequalities 

< x" = ~ (2a 2 + b 2 - xx/2 + 2/3-7) < 1 
< 2 Ret' = ~(ai + 2bx - xi/2 - (j. - 2v) < 1 



■Ea - 



5S 2 



(86) 



(87) 



•1 < a+ g(2ai - 61 -xi - 2fj, + v) < 1 



(89) 



are satisfied. Otherwise the summand is taken to be zero. 

Again, in a similar fashion to the previous section, we note that if the above approximation is to become O(l), 
we must set 2a± — b\ — x\ — 2/i + v = and a 2 — 2b 2 + x\ + (3 + 27 = 0. Consequently, under the assumption 
< at, a 2 , &i, &2 < 1) the above three inequalities will now require us to set fj, = v = a = 8 = 7 = and x\ = [2ai\ . 
Furthermore, we are assuming ai ^ 1/2 and b 2 ^ 1/2 by ignoring cases of equality in (87-89J). 

Thus, by discarding all exponentially small terms, we arrive at the semi-classical propagator 



(b\B e ,s 



~(l - i|2oi - h - L 2a iJ I ) ,>x i> 



- ^{(2ai - 61 - L2aiJ) 2 + (03 - 2b 2 + [2 ai \Y 



i(3aia 2 + 36i6 2 + 4ai& 2 - Aa 2 b x ) - 2i[2a x \ (a x + 2b x + 2a 2 + b 2 - [2aiJ/2)| 



• exp 



■kD 



{(2ar - h - L2 fll J) 2 + (a 2 - 2b 2 + [2 ai \Y 



i(3a!a 2 + 3&1&2 + 4a x 6 2 - 4a 2 6x) - 2«L2a!j (a a + 2b x + 2a 2 + 6 2 - [2a a J /2) V 



+ o(l) (90) 



+ o(l). (91) 



Both forms, (|90j) and (|9lj), are equally valid since their difference is exponentially small (although the first (|90|) may 
prove to be more accurate). Comparing (^lj) to (f50|), we see that for < 9 < 1, the semi-classical propagator also 
takes the Van Vleck form ([5l|). Furthermore, the classical baker's map will be recovered in the limit D — > 00. 



Case = 1. 



After the rousing success of the previous calculations, it is tempting to conclude that the classical baker's map will 
always be restored in the limit D — > 00. Unfortunately, when 6 — 1 (0 = 0), certain assumptions made previously 
will prove incorrect. In particular, in Eq. (^2|), we have assumed that terms will dominate; however this clearly 
cannot now be the case. In this section we show that such differences cripple any hope that the classical baker's map 
will be recovered for all possible classical limits. 



When 9 = 1 the number of momentum bits r = — s > remains constant. Using (|16|) and (17), with n = N — r 
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and R = 2 r , our matrix element is 



oo 1 fl-l 2R-1 D/{2R)-X 

(b\B x ,- r \a) = * E E E E E exp -^(H 2 + |6| 2 -a 2 -6* 2 )+z.( M -,) 
- TTD(Rq m + x x /2 - a + fi + (k + 1/2 - R/2)/D) 2 - nD(2Rq m - b* + v + (l + 1/2 - i?)/L> 



+ ^{ (i + 1/2) (i + 1/2) + iter (J + 1/2) -2(j + 1/2) (fe + 1/2) } 



(92) 



D 



R~l 2R 



-1 1/2 



E E E E /* ex p -^(H 2 + H 2 - a 2 -^)+^-,- K ) 

t ™ n a i» n I n « 



/x, f ,K— — oo a; i—0 /— 

- 7rL>(i + a;i/2 -a + ^i) 2 -7rD(2i- fe* +^) 2 - -^(fc + 1/2 - i?/2) 2 - ^(l + 1/2 - i?) 2 
2«7rD 

+ — — *k - 27r(t + xi/2 - a + /i) (k + 1/2 - R/2) - 2ix{2t - b* + v) (l + 1/2 - R) 

+ J[{ti+ V2) (i + 1/2) + iZsi (I + 1/2) - 2{j + 1/2) (k + 1/2) } (93) 

where we have converted the sum over m to an integral over t using the same technique in the previous sections. 
Now defining the integral 

1/2 

h= I dt exp \-irDf{t)-2n(t-A)(k + l/2-R/2) - 2n(2t - B)(l + 1/2 - R)] (94) 



where 



one finds the saddle point 
and hence, the approximation 

oo+2 Im t' 



f(t) = (t-A) +(2t-B) -2Uk/R 
A = a — x\/2 — /i 
B = b* - v 



A + 2B in 



5R 



h 



dt exp - 7rDf(t) - 2%(t - A)(k + 1/2 - R/2) - 2n(2t - B) (l + 1/2 - R) 



-oo+i Im t' 

l 

—== exp 

2vr 



nD 



, 2 2iirD 



■nD 



5i? 2 



2z7T 7T 9 

— |2A - B) (2k - I + 1/2) — (/c + 2/ + 3/2) k H — — (fc + 2Z + 3/2 — 5R/2) 



(95) 
(96) 
(97) 



(98) 



(99) 



(100) 



provided that < Ret' < 1/2, and otherwise zero. 

Apart from the last term in the exponent, Eq. (100) is simply the well-known formula for a saddle-point approx- 
imation found in any standard text. We will now drop this 1/D term, along with all others in (193), to obtain the 



15 



approximation 

(b\B lt _ r \a) 



V5i? 2 



oo 1 R-l 2R-1 

E E E E - P 

H,u,k=— oo X!=0 j,k=0 1=0 



ttD r, 



b\ — x\ — 2/i + lA 



if 



+ (a 2 — 26 2 + K /-R) 2 | + wr(/-« — v) + ii^R)(a\a 2 — 6162) 

— |(2ai - 61 - xi - 2/i+i/) (2a 2 + 6 2 ) - (ai + 2&i - x x /2 - /i - 2i/)k/.RJ 

27T 277T 

+ — (2A- B) (2k- 1 + 1/2) - -— (fc + 2Z + 3/2)k 

D oil, 
ITT ( ^ 

+ { (i + 1/2) (/ + 1/2) + ifcci (Z + 1/2) - 2(i + 1/2) (fc + 1/2) } 



< 2Rei' = + 2b\ - si/2 - fx- 2v\ I 



(101) 



(102) 



and otherwise zero. 

As in the previous cases, under the assumption < 01, 61 < 1, we may use ( |l02j ) to set fx = v = and = [2aiJ . 
Thus, we arrive at the following semi-classical approximation for our propagator 



(b\B x ^ r \a) 



1 



oo R-l 2R-1 

E E E ex p 

e=— oo j.k=0 1=0 



h- L2aiJ) 2 + (a 2 -2b 2 + n/R) 2 



+ i(3aia 2 + 36i6 2 + 4ai6 2 - 4a 2 6i) - 2in(a 1 + 2b x - [2a 1 \/2)/R - 2z[2aiJ (2a 2 + 62) j 

2ir 2iir 
+ —(2a- b* - [2aiJ)(2fc- 1 + 1/2) (fc + 2Z + 3/2) K 

+ ^{(j + 1/2) (l + 1/2) + R[2 ai \ (I + 1/2) -2(j + 1/2) (fc + 1/2) } 



0(1) 



(103) 



Note that the summation index k remains unconstrained. Consequently, additional probabilistic 'humps' emerge at 
locations other than those specified by the classical baker's map. In fact, in the region < 61,62 < 1, there will be 
2R humps at the positions (61, 62) = (2ai — \ 2a\\ , (da + n/R)/2) where —a 2 R < k < 2R— a 2 R. 
Consider the simplest case r — 0. Our semi-classical propagator is then 

exp - ZL-j(2ai -&i - L2aiJ) 2 + (a2 - 26 2 + K) 2 +i(3aia 2 + 36i6 2 +4ai62-4a 2 6i) 
k=o L 

- 2i«(o 1 + 26i - L2ajJ /2) - 2i [2a x \ (2a 2 + 6 2 ) | + wr( L2aiJ - k) 

o(l) 



cos 



— (2a- 6* - L2aiJ) + -(L2a x J + 1/2) - — n 

o ZD 



(104) 



defining two humps: one at a position specified by the classical baker's map, (bt, b 2 ) = (2a\ — [2ai\ , (a 2 + |_2aiJ)/2), 
with an asymptotic size of 



\(b(a)\B 1<0 \a) 



(a 2 - 1/2) 



and another at (61,62) = (2ai — [2aiJ , (02 + 1 — |_2aiJ)/2) with the size 



|<6(a)|Si j0 |a) 



■(as -1/2) 



o(l) 



o(l) 



(105) 



(106) 



One interpretation of these equations could be that a stochastic mapping is implied in the classical limit: a point at 
(01, a 2 ) has the probability cos 2 [7r(a2 — 1/2)/2] of obeying the classical baker's map, and probability sin 2 [ir(a 2 — 1/2)/2] 
of ending up at (2ai — [2aiJ , (02 + 1 — |_2 a iJ )/2). Notice that there is now a smooth transition of probabilities as one 
crosses the lines a 2 = 0, 1. 



16 





Consider the size of our probabilistic humps in the general case, 
with k fixed, then 



If we set (61, 6 2 ) = (2oi - L2aiJ , (a 2 + k/R)/2), 



(b(a)\B lt _ r \a) = 



-exp 



iirD 
~2R 



where 



R-l 



[2a\n + R[2ai\a 2 — [2aiJ k 



[■jTRa 2 ] 



tf„(a)+o(l) 



=0 sin 



7r(a 2 -(i + l/2)/^ 









): 


sin 











a 2 - L2aiJ + k/R - (j + 1/2)/ R 



(107) 



(108) 



The probability associated with each hump is given by $>l, with E«=o *k = ^ and is plotted for r = 2 in Fig. 
|(a). The curve with the largest probabilities is that associated with the 'correct' hump prescribed by the classical 
baker's map («' = |_2aiJ-R). It takes maximum values of unity whenever a 2 = (m + 1/2) / R (0 < m < R— 1), and as 
the number of momentum bits (r = log 2 R) becomes large, ^> 2 K , — > 1 for all < a 2 < 1. Hence, as expected, we are 
left with a single hump located at the correct position. This agrees with the < 9 < 1 case. In figures ||(c) and (d), 
we have drawn the function |(6|f?|a)| 2 and its semi-classical approximation (103), respectively, when a = 3/4 + «/20, 
r = 2 and N = 8. For such a large dimension, D = 2 8 , our semi-classical approximation becomes almost identical 
to the exact matrix element. One may also view our semi-classical propagator as an approximation to the Husimi 
function for the mapped state B\a). 
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IV. CONCLUSION 



In this paper, we have derived the semi-classical form of the Schack-Caves quantum baker's maps, taking into 
consideration the different possible ways of obtaining the classical limit. We have shown that whenever the number of 
momentum bits becomes infinite in the limit Ti — > 0, the quantum propagator in the coherent-state basis takes the Van 
Vleck form (|5l|). Therefore, we conclude that the classical baker's transformation is restored for all such cases. In the 
case where the number of momentum bits is held constant, the classical limit is not that of the baker's transformation, 
but a stochastic variant. It may be possible one day to explore this discrepancy experimentally with the help of a 
quantum computer [Q, |). As a final note, we remark that although our semi-classical formula (|m]) espouses a certain 
familiarity, we should not be presumptuous about its generality. One should, in general, expect extra phases in the 
exponent (see Baranger et al j41j). This has already been found in a spherical geometry for the case of the kicked top 
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